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This  report  describes  the  second  year's  work  on  this  research  contract 
which  has  the  aim  of  developing  and  utilizing  measurement  and  data  analysis 
techniques  to  investigate  the  structure  of  turbulent  reacting  flows  with 
well  defined  initial  and  boundary  conditions.  Experiments  are  in  progress 
in  three  fields:  free  jet  diffusion  flames,  flames  impinging  on  flat 
plates  and  turbulent  fuel  sprays.  The  major  initial  effort  has  concentrated 
on  the  jet  flame  and  new  and  significant  information  has  been  derived  on 
coherent  eddy  structures  in  the  initial  regions  of  flames.  Experimental 
investigations  of  aspects  of  diffusion  flame  structure  have  also  received 
financial  support  from  the  Office  of  Naval  Research  (Project  SQUID) , 
Subcontract  No.  8960-30.  This  research  has  also  benefited  from  the 
assistance  of  personnel  in  the  Combustion  Aerodynamics  Research  Laboratory 
funded  by  US  Army  Grant  DA-ER0-79-0031  for  the  study  of  impinging  jet 
flames  (from  4/1/79)  and  NASA  Research  Grant  -  No.  NSG-7517  for  the  study 
of  vaporizing  fuel  sprays.  In  addition  use  has  been  made  of  major  equipment 
items  and  personnel  funded  by  the  Science  Research  Council  (U.K.)  and  the 


INTRODUCTION 


The  objective  of  this  research  is  to  obtain  increased  understanding  of 
turbulent  combustion  in  flows  related  to  propulsion  systems.  Measurement  and 
data  analysis  techniques,  and  the  boundary  conditions  of  various  flames,  are 
designed  to  provide  an  improved  conceptual  framework  and  fundamental 
experimental  data,  to  aid  the  development  and  checking  of  modelling  techniques 
for  practical  combustion  systems.  The  specific  approach  of  the  study  is 
to  quantify  the  roles  of  large  eddies  (coherent  structures)  in  reacting 
flows.  The  experiments  and  measurement  techniques  which  have  been  developed 
to  attain  this  aim  were  described  in  the  previous  Interim  Scientific  Report. 
Measurements  during  the  past  year  have  been  concentrated  on  the  initial 
region  of  a  gaseous  jet  flame  with  separate  variation  in  Reynolds  number 
and  equivalence  ratio.  These  measurements  have  produced  interesting  and 
detailed  data  on  the  structure  of  the  initial  regions  of  jet  diffusion 
flames,  as  a  function  of  systematic  variation  of  initial  burner  nozzle 
conditions.  This  data  is  already  recognized  to  be  of  considerable  poten¬ 
tial  use  to  combustion  modelling  groups. 

Two  associated  experiments,  investigating  the  structures  of  turbulent 
liquid  fuel  sprays  and  impinging  flames  have  been  initiated  during  the 
past  year.  These  experiments  will  also  concentrate  on  the  roles  of  eddy 
structure  in  controlling  reaction  in  these  types  of  flows.  The  present 
report  concentrates  on  the  extensive  fundamental  data  which  has  been 
derived  on  the  structures  of  the  jet  flames. 
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JET  FLAME  EXPERIMENT  -  PHILOSOPHY  AND  DESIGN 

Figure  1  shows  the  arrangement  and  notation  for  the  jet  flame 
experiment.  The  axisymmetric  flame  is  produced  by  the  flow  of  fuel  gas 
through  a  25.4mm  diameter  nozzle  into  a  400mm  x  400mm,  low  velocity  (0.73-lm/s), 
secondary  air  flow. 

The  experiments  are  designed  to  derive  time  dependent  information 
from  measurements  of  fluctuating  quantities  in  diffusion  flames.  Fundamental 
information  and  insight  are  required  on  the  'physics’  of  reacting  flows 
and  this  information  must  be  sufficiently  extensive  and  accurate  to  be  of 
use  to  those  involved  in  attempting  to  model  and  predict  turbulent 
combustion.  A  special  effort  is  made  in  these  experiments  to  control 
carefully,  vary  systematically  and  measure  the  flame  initial  and 
boundary  conditions.  By  this  means  detailed  measurements  made  in  the 
flame  can  be  compared  directly  with  modelling  predictions  based  on  these 
boundary  conditions.  The  burner  nozzle  is  contoured  to  give  a  flat 
velocity  profile  at  the  exit  with  laminar  internal  and  external  boundary 
layers.  The  initial  turbulence  levels  of  both  the  primary  (fuel  jet)  and 
secondary  flows  are  less  than  1%.  Ignition  and  stabilization  take  place 
in  a  laminar  region  near  the  nozzle  lip.  By  progressively  making  measurements 
from  the  initial  conditions  at  the  nozzle  exit,  the  natural  development 
of  the  flame  is  followed,  from  the  laminar  conditions  at  the  nozzle  exit, 
through  a  transition  region,  to  the  fully  turbulent  flow  downstream.  This 
experimental  forward  marching  procedure  will  lead  to  fundamental  under¬ 
standing  of  the  physical  phenomena.  Diagnostic  techniques  have  been, 
and  continue  to  be,  developed  for  making  measurements  with  high  spatial 
and  temporal  resolution. 

The  combination  of  high  speed  visualization  and  measured  data  is 
already  beginning  to  show  important  deterministic  features  of  the  flows. 

Thus,  features  such  as  burning  interface  layers,  modes  of  stability  of 


Interfaces,  reacting  vortices,  coherent  luge  eddies  etc.,  will  be 
investigated  by  suitable  measurement  and  data  analysis  techniques. 


mm 
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(v)  High  Speed  Cine  Photography 

The  high  speed  cine  camera  provides  detailed  qualitative  pictures 
of  the  flame  structure.  Frame  by  frame  analysis  of  the  films  provides 
quantitative  information  on  the  dimensions,  velocities  and  passing 
frequencies  of  'events'  in  different  regions  of  flames.  The  photographic 
information  also  provides  the  essential  framework  for  the  subsequent 
interpretation  of  fixed  point  measurements  in  terms  of  flame  structure, 
where  such  structure  has  recognisable  organisation.  Most  of  the  High 
Speed  Films  have  been  taken  using  a  colour  schlieren  system  with  600mm 
diameter  mirrors. 

EXPERIMENTAL  RESULTS 

The  bulk  of  initial  measurements  have  been  made  in  a  propane  flame. 

Jet  Flame  I,  under  the  following  initial  conditions  at  the  nozzle  exit: 

4 

jet  velocity,  =  6.4m/s;  Reynolds  Number,  Re  =  10  ;  propane/air 
equivalence  ratio,  $  =  10.4.  Measurements  have  been  made  in  the  first 
20D  of  flame.  Measurements  have  also  been  made  with  variation  in  Re  and  $. 

The  1978  scientific  report^  described  the  results  of  direct  observa¬ 
tions  of  the  flame  based  on  flame  luminosity.  The  occurrence  of  vortex 

motions,  distorting  the  flame,  and  the  coalescence  of  these  vortices 

\ 

beyond  X  =  16D  was  noted.  However  the  more  recent  schlieren  films  and 
extensive  point  measurements  have  now  revealed  a  flow  structure  which  is 
considerably  more  complex  them  was  first  assumed.  This  complexity  is 
caused  by  the  simultaneous  occurrence  of  several  modes  of  instability 
near  the  nozzle  so  that  more  than  one  type  of  'coherent  structure'  can 
be  observed  at  any  one  axial  position.  It  has  been  found  that  the  Kelvin- 
Helmholtz  instabilities,  which  were  at  first  assumed  to  dominate  the 
transitional  flame  are  accompanied  by,  equally  important,  combustion 
instabilities  which  are  mainly  influenced  by  the  physico-chemical  properties 
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of  the  fuel  gas  mixture,  rather  than  the  fluid  dynamics. 

Flow  Visualizations 

Colour  schlieren  cine  films  have  been  made  of  a  series  of  flames 

with,  in  the  first  instance,  systematic  variation  in  equivalence  ratio, 

$,  with  Reynolds  number.  Re,  fixed.  Subsequently  $  was  systematically 

varied  while  Re  was  maintained  constant.  This  allows  a  separate  study 

of  the  effects  of  fluid  mechanical  and  combustion  mechanisms  in  the 

transitional  flames.  The  term  'transition*  refers  to  the  transitional 

flow  between  the  laminar  region  near  the  nozzle,  to  fully  turbulent  flow 

further  downstream.  The  Reynolds  number  (Re  =  U  D/v  )  has  been  varied 

J  AIR 

NTP 

3  4 

between  3  x  10  and  1.5  x  10  .  The  equivalence  ratio,  4>,  has  been  varied 

between  “  ('pure'  diffusion  flame)  and  1.5  (i.e.  approaching  premixed 

conditions) .  In  all  of  the  flames  studied  vortex-like  structures  could 

be  detected  in,  at  least,  the  first  20D  of  flow.  Visually,  these  structures 

appear  to  be  similar  to  those  reported  earlier  '  '  by  the  authors. 

Significant  changes  can  be  observed  when  both  $  and  Re  are  varied.  Most 

of  the  data  that  has  been  analysed  to  date  is  for  ’Jet  Flame  I’,  with 
4 

Re  =  10  and  <p  *=  10.4.  A  general  description  of  visualisations  of  this 
flame  is  given  below.  Quantitative  data  from  the  films  is  compared  with 
point  measurements  later  in  this  report. 

The  sketches  shown  in  Figures  2,  3  and  4  have  been  derived  from 
frame  by  frame  observations  of  the  films  and  describe  the  general  appearance 
of  the  first  15D  of  Jet  Flame  1,  divided  for  convenience  into  5D  sections. 

The  most  obvious  feature  seen  on  the  films  is  a  'double*  structure  in  the 
first  20D,  with  inner  fast  moving  eddies  (Convection  Velocity  =  6m/s)  and 
outer,  large,  slow  moving  eddies  (Convection  Velocity  <  2m/s) .  The  inner 
eddies  have  the  appearance  of  hot/cold  interfaces  developing  wave  instabilities. 
These  interface-waves  coalesce  and  interact:  wave  amplitudes,  wave  lengths 
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and  spacings  between  waves  all  Increase  with  increasing  distance  downstream. 

The  outer  vortex-like  eddies,  which  are  also  indicated  by  hot/cold  interfaces, 

develop  much  more  slowly.  There  appears  to  be  no  observable  correlation 

or  interaction  between  the  outer  and  inner  vortices  up  to  about  x  =  lOD. 

However,  beyond  this  position,  the  outer  eddies  appear  to  encroach 

increasingly  into  the  centre  of  the  jet  and  there  is  an  obvious  complex 

interaction  between  these  two  instability  modes.  The  outer  hot  flow 

region  has  been  termed  the  'pre-heat  zone'  (Figs.  2,  3  and  4)  and  the 

inner  interface  has  been  termed  the  'reacting  interface'.  This  pre-heat 

zone  can  be  clearly  seen  in  the  classical  photograph  of  Wohl,  Gazley  and 
4 

Kapp.  More  precise  information  and  a  more  complete  description  of  the 
main  reaction  zones  awaits  the  acquisition  of  more  extensive  data,  using 
conditional  sampling  techniques,  from  measurements  using  the  ionization 
and  thermocouple  probes. 

With  increasing  distance  downstream  there  is  an  increasing  growth  of 
'three-dimensionality'.  This  is  observed,  first  by  verticil  'streaks', 
and  later  by  'cellular'  interface  structures,  reminiscent  of  schlieren 
photographs  of  turbulent  premixed  flames  (e.g.  Lefebvre).  Photographs  taken 
by  direct  photography,  using  flame  luminescence,  show  that  the  initial 
streamwise  streaks  cure  similar  to  those  observed  in  cold  jets**  in  which 
vortex  rings  develop  azimuthal  waves  as  one  of  the  stages  in  the  transition 
process.  Films  of  flames  with  varying  <J>  indicate  that  the  local  thickness 
of  the  'preheat  zone'  increases  as  $  decreases;  as  the  jet  mixture 
ratio  is  made  leaner,  the  stoichiometric  mixture  'line'  moves  further  in 
towards  the  jet  centre.  Varying  $  also  changes  the  frequencies  of  the 
inner  vortices;  this  important  effect  is  discussed  in  the  section  des¬ 
cribing  the  laser-schlieren  measurements.  Increasing  Re  resulted  in  the 
more  rapid  onset  of  three-dimensionality,  and  thus,  presumably,  a  reduction 
in  the  transition  length  of  the  flame. 
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Interpretation  of  the  schlieren  films  was  hampered  by  the  inability 
to  interpret  line-of-sight  integrated  data  unambiguously  in  terms  of 
point  information.  At  first  sight,  the  films  give  the  impression  of 
providing  a  cross-sectional  view  of  the  flame  in  a  radial  plane,  perpendicu¬ 
lar  to  the  optical  axis.  The  initial  laminar  flame  has  a  thin  cylindrical 
shape  so  that  most  of  the  observed  information  indicates  the  flame 
structure.  As  instability  develops,  radiation  from  the  'front'  and  'back' 
parts  of  the  flame  is  recorded  on  the  film.  Further  insight  into  this 
interpretation  problem  will  be  gained  by  making  use  of  the  probe  measure¬ 
ments  and,  in  particular,  by  filming  the  flow  at  the  same  time  as  measure¬ 
ments  are  made. 

Velocity  Measurements 

Cold  Jet  I,  the  air  jet  with  the  same  Re  value  as  Jet  Flame  1,  was 

studied  in  detail  by  using  a  hot  wire  anemometer.  The  nozzle  boundary 

2  -1 
layer  is  laminar  with  a  'shear  thickness'  U(6U/6r)  =  2.36mm.  The 

J  max 

nozzle  exit  turbulence  level  is  0.4%,  and  these  small  fluctuations  occur 
mainly  at  a  frequency  of  30  Hz,  which  is  assumed  to  be  the  resonant  fre¬ 
quency  of  the  jet  settling  chamber.  The  initial  instability  frequency  of 
the  mixing  layer  is  215  Hz  which  is  within  5%  of  that  predicted  by  the 

7 

theory  of  Michalke.  Hot  wire  frequency  spectra  in  Cold  Jet  I  showed  the 
appearance  of  subharmonics,  and  thus  a  vortex  ring  coalescing  region, 
between  *  =  2D  and  x  »  4D.  The  frequency  spectra  approach  the  turbulent  form 
beyond  x  =  4D,  indicating  that  this  is  the  transition  length  of  the  cold  jet. 

Figure  5  shows  LDA  measurements  of  mean  velocity  in  Cold  Jet  I.  The 
data  is  in  excellent  agreement  with  hot  wire  measurements,  and  also  with 
previously  reported  measurements.  This  indicates  the  good  general  accuracy 
of  the  LDA  system.  This  data  was  obtained  without  frequency  shifting  so 
that  inaccuracies  are  possible  (particularly  for  fluctuating  velocity  data) 
where  the  mean  velocity  falls  below  2m/s.  The  low  velocity  regions  of  the 
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flows  will  be  measured  in  detail  by  using  a  Bragg  Cell  (shifted)  LDA. 

The  potential  core  of  Cold  Jet  1  ends  at  I  =  4D. 

Figure  6  shows  mean  velocity  profiles  in  the  first  20D  of  Jet  Flame 
1.  These  results  are  in  dramatic  contrast  to  the  data  shown  in  Figure  5, 
for  the  non-burning  flow  at  the  same  Reynolds  number.  The  potential  core 
in  the  flame,  as  indicated  by  the  central  uniform  velocity  region,  now 
apparently  extends  beyond  x  =  20D.  The  centre  line  velocity  changes 
little  in  magnitude  in  the  first  20D  length  of  flame.  Mean  velocity 
gradients  are  generally  higher  in  the  flame,  than  at  the  same  axial 
position  in  the  cold  jet.  However  the  total  width  of  the  flows  do  not 
differ  greatly  between  the  flame  and  cold  cases,  for  any  particular  x 
value.  Examination  of  individual  velocity  profiles  in  the  flame  reveals, 
at  some  stations,  the  occurrence  of  ‘humps'  in  the  profiles,  i.e.  more 
than  one  inflexion  point  is  found  compared  with  the  classical  single 
inflexion  point,  error  function  or  Gaussian,  shapes  in  the  cold  jet.  These 
and  other  differences  between  the  cold  and  burning  flows  need  to  be  explained 
in  terms  of  the  various  phenomena  introduced  by  chemical  reaction  and 
heat  release.  Attention  is  being  focussed  initially  on  the  effects  of 
combustion  on  the  relatively  orderly  structure  of  the  transitional  flows; 
in  the  first  place  to  obtain  a  qualitative  understanding  of  the  various 


processes,  and  secondly  to  derive  quantitative  data  on  these  individual 
processes  for  eventual  comparison  with  predictions  with  the  aim  of 
developing  more  accurate  computer  modelling  techniques. 

— r  h 

Figures  7  and  8  show  profiles  of  fluctuating  velocity,  (U)  /U  , 

J 

measured  by  LDA  in  Cold  Jet  1  and  Jet  Flame  1  respectively.  A  special 
experimental  problem  arises  in  measuring  turbulence  intensity  in  the 
vicinity  of  peaks  in  the  region  very  near  the  nozzle  (X  =  O.Q4D) .  In 


this  region,  an  anomaly  is  introduced  by  the  relatively  large  LDA  measure- 
volume 


ment^when  compared  with  the  shear  layer  thickness. 


LDA  measurements  do 
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not  agree  with  hot  wire  data  which  show  turbulence  levels  less  than  1%  across 
the  complete  flow  in  this  region.  In  common  with  most  published  investiga¬ 
tions  using  LDA,  the  LDA  turbulence  levels  tended  to  be  higher  than  those 
measured  by  hot  wires  by  a  factor  of  U  /lOO,  approximately.  Thus  a  1.5% 
turbulence  level  measured  by  LDA  corresponded  to  a  0.5%  level  measured 

by  hot  wires.  This  overestimation  of  turbulence  level  by  LDA  is  a  result 

8 

of  an  accumulation  of  errors  and  biasing  inherent  in  all  LDA  systems. 

Further  attempts  are  being  made  to  improve  the  signal  processing  and 
validation  techniques  so  as  to  provide  more  accurate  LDA  turbulence 
measurements.  At  this  stage  we  are  confident  that  data  presented  in 
this  report  is  at  least  comparable  in  accuracy  with  the  best  of  previously 
published  LDA  measurements  in  flames. 

The  turbulence  distributions  in  Figure  7  are  in  good  agreement  with 
those  reported  previously  in  cold  round  jets,  using  both  hot  wires  and 
LDA.  Peak  turbulence  levels  occur  near  the  end  of  the  potential  core 
at  x  -  4D  followed  by  a  decay,  approaching  the  slope  x  1  for  exact 
similarity.  The  flame  core.  Figure  8,  is  quite  different,  with  a  much 
slower  increase  in  turbulence  levels  with  increasing  distance  downstream. 
However  the  maximum  levels,  at  x  =  20D,  are  higher  than  those  found  in 
the  cold  jet,  i.e.  20%  as  opposed  to  15%.  Furthermore  there  is  every 
likelihood  of  further  increases  beyond  x  =  20D,  a  region  in  which  measure¬ 
ments  have  not  yet  been  made.  The  peak  turbulence  levels  in  the  flame 
are  found  nearer  the  outer  low  velocity  stream  and  there  cure  relatively 
small  turbulence  levels  along  the  centre  line  of  the  mixing  layer  at  the 
nozzle  lip  radius.  Double  peaks  in  the  turbulence  distributions  are 
seen  up  to  x/D  =  6  in  Figure  8.  These  could  be  connected  with  the 
interactions  of  vortices  or  waves  in  the  transitional  flame.  These  results 
show  that  there  is  a  significant  influence  of  combustion  on  the  turbulence 


field  in  the  initial  region  of  the  flame.  Further  work  is  required  in 


L3 


order  to  explain  these  findings  and  develop  models  which  are  consistent 
with  these  phenomena. 

In  two-stream  mixing  flows,  such  as  those  investigated  here,  it  is 
necessary  to  provide  uniform  seeding  in  both  streams  so  as  to  avoid 
biasing  effects  in  LDA  measurements.  The  influence  of  changing  seeding 
particle  number  density  in  each  stream,  nas  not  been  reported  in  any 
detail  in  the  literature.  A  series  of  tests  were  carried  out  in  which 
the  relative  seeding  densities  of  the  primary  and  secondary  flows  were 
progressively  varied.  Very  significant  changes  could  be  produced  in 
measured  mean  and  r.m.s.  velocities,  depending  upon  the  relative  seeding 
levels,  and  the  position  in  the  flow.  This  indicated  the  need  for  careful 
equalisation  of  seeding  levels  before  each  set  of  experiments.  Separate 
from  the  problem  of  making  accurate  measurements  of  time  averaged  quantities 
controlled  changing  of  seeding  levels  can  be  used  to  derive  conditionally 
sampled  velocity  data  which  cam  be  interpreted  in  terms  of  intermittency 
and  fuel/air  mixing.  Figure  9  showc  am  example,  in  Jet  Flame  1,  in  which 
mean  velocity  profiles  are  measured  with,  (i)  primary  flow  seeding  only, 

(ii)  secondary  flow  seeding  only  and  (iii)  both  flows  seeded.  In  general 
the  'primary  seeding'  velocity  is  higher  than  the  local  mean  velocity, 
and  the  'secondary  seeding'  velocity  is  lower.  This  agrees  with  what 
could  be  expected  from  intermittency  sampled  data  obtained  in  mixing  layers 
in  the  past,  by  using  hot  wires.  It  is  interesting  to  note  that  the 
primary  flow  (including  the  products  resulting  from  burning  of  this 
primary  flow)  penetrates  only  to  the  centre  of  the  mixing  layer.  The 
secondary  air  penetrates  the  full  width  of  the  mixing  layer  up  to  the 
edge  of  the  potential  core.  Mixing,  and  thus  reaction,  can  only  occur 
in  the  region  of  overlap  of  the  two  streams,  i.e.  approximately  7mm  <r<  18mm 
at  x  =  6D.  It  will  be  shown  that  this  corresponds  to  the  region  of  high 


ionization  levels. 


Ionization  Levels 


Measurements  of  ionization  levels  have  been  made  in  Jet  Flame  1, 
including  mean  values,  r.m.s.,  p.d.f.'s  and  frequency  spectra.  Problems 
of  calibration  drift,  condensation  and  sooting  of  the  probe  have  been 
encountered;  these  create  particular  difficulties  with  regard  to  fluctuating 
ionization  density  measurements.  These  problems  have  now  been  largely 
overcome  in  the  new  designs  of  probe.  Distributions  of  mean  ionization 
level  are  shown  in  Figure  10.  It  is  seen  that  the  region  of  significant 
ion  density  is  relatively  narrow,  compared  with  the  total  mixing  layer 
width  indicated  by  the  mean  velocity  distributions.  High  ionization 
levels  indicate  high  rates  of  reaction  at  elevated  temperatures.  Compari¬ 
son  of  Figures  9  and  lO  shows  that  the  regions  of  overlap  of  the  primary 
and  secondary  flow  seeding  particles  (where  mixing  takes  place)  corres¬ 
pond  with  regions  of  high  ionization  levels.  Peak  signals  are  found  near 
x  =  6D.  Interpretation  of  signal  magnitude  is  not  straightforward 
because  of  the  complex  dependence  on  probe  geometry,  flow  velocity  and 
reacting  layer  width,  as  well  as  ion  density.  The  sensitivity  of  the 
ionization  probe,  to  the  occurrence  of  reaction,  makes  it  a  valuable  tool 
for  studying  the  transitional  flame,  where  convoluted,  reacting  interface 
layers  are  clearly  visible. 

LDA  mean  velocity  and  ionization  probe  measurements  at  x  =  8D  in 
Jet  Flame  1  are  compared  in  Figure  11.  It  is  interesting  to  note  that 
there  is  a  'hump'  in  the  mean  velocity  profile  near  the  peak  ionization 
level  position.  This  can  be  explained  as  being  due  to  flow  acceleration 
accompanying  reaction  and  heat  release. 

Temperature  Measurements 

Temperature  measurements  have  concentrated  initially  on  a  detailed 
mapping  of  the  mean  temperature  field  of  Jet  Flame  1  using  fine  wire 
thermocouples  with  digital  processing  and  radiation  corrections. 


Measurements  of  fluctuating  temperature  are  currently  in  progress. 

Figure  12  shows  mean  temperature  profiles  in  Jet  Flame  1.  The  peak  mean 
temperature  is  found  near  x  =  3D.  The  temperature  peak  moves  progressively 
away  from  the  centre  of  the  jet  flame  and,  beyond  x  =  3D,  decays  in  value 
with  increasing  distance  downstream.  The  peak  temperature  is  always 
further  from  the  jet  centre  than  the  peak  ionization  level  (c.f. 

Figs.  10  and  12).  The  separation  between  the  two  peaks  increases  with 
increasing  distance  downstream.  The  reason  for  this  difference  is  not 
certain  but  physical  models  which  should  explain  the  phenomena  are  being 
formulated  and  will  be  checked  by  further  experiments.  At  each  station, 
the  total  local  mixing  layer  width,  indicated  by  the  mean  temperature 
data,  lies  within  lO%  of  the  width  indicated  by  the  velocity  profiles. 
Laser-Schlieren  Data 

Laser-schlieren  spectra  have  been  measured  at  different  longitudinal 
positions  in  Jet  Flame  1,  with  the  laser  beam  passing  through  the  flame, 
tangential  to  a  cylindrical  projection  of  the  nozzle  lip  (r  =  0.5D) . 

These  spectra  represent  integrated  data  for  the  complete  length  of  the 
beam,  and  they  thus  contain  contributions  from  values  of  r  greater  than 
0.5D.  The  transformation  of  this  integrated  line  of  sight  data,  into 
'point'  measurements  is  being  investigated  by  using  an  Abel  transformation. 

Figure  13  shows  spectra  measured  at  four  axial  stations.  Distinct 

peaks  are  found,  corresponding  to  the  regular  wave  and  vortex  structures 

seen  in  the  schlieren  films.  Several  distinct  peaks  are  found  at 

each  axial  station;  in  general,  there  are  no  clear  harmonic  relationships 

between  the  frequencies  at  which  these  peaks  occur.  Peak  amplitudes  of 

the  distinct  frequencies  rise  to  maximum  values  at  x/D  -  4  and  then  decline 

with  distance  downstream.  This  is  in  some  accord  with  the  appearance  of 

subharmonics  of  the  initial  Kelvin-Helmholtz  instability  in  spectra  in 

2 

the  transitional  cold  jet.  However  there  are  important  and  significant 


differences  between  the  cold  and  flame  cases;  this  is  made  clear  by 
the  results  of  a  series  of  systematic  experiments  in  which  Re  and  $  were 
varied. 

Figure  14  shows  spectra  obtained  when  Re  was  fixed  and  the  equivalence 
ratio  $  was  varied.  It  can  be  seen  that  there  is  a  gradual  decrease  in 
the  most  energetic  frequency  as  $  was  increased  by  increasing  the  propor¬ 
tion  of  fuel  in  the  propane/air  mixture,  while  maintaining  a  constant 
exit  velocity.  Similar  results  were  obtained  at  higher  Reynolds  numbers, 
but  the  most  energetic  peaks  then  occurred  at  higher  frequencies.  Figure  15 
shows  spectra  obtained  by  maintaining  a  fixed  equivalence  ratio  and  varying 
the  Reynolds  number  (i.e.  varying  the  primary  flow  rate).  It  is  seen 
that,  for  all  Re,  there  are  a  series  of  distinct  peaks  in  the  spectra, 
but  the  frequency  values  of  these  peaks  are  independent  of  Reynolds 
number,  it  is  deduced  that  the  inner  vortices,  which  distort  the  initially 
laminar  reacting  interface  layer  (Fig.  2) ,  cannot  be  simply  explained  as 
being  due  to  Kelvin-Helmholtz  Instabilities,  as  generally  found  in  cold 
jets,  even  if  account  is  taken  of  heat  release.  Rather,  it  appears  that 
combustion-driven  instabilities  are  generated  which  result  in  wave  and 
vortex  structures  similar  in  appearance  to  the  Kelvin-Helmholtz  vortices. 

Thus  there  are  a  series  of  discrete  frequencies  at  which  the  inner 
instability  cm  occur  and  the  values  of  these  frequencies  depend  on  the 
physico-chemical  properties  of  the  fuel  gas.  The  flow  rate  (i.e.  Re) 
determines  the  frequency,  or  batch  of  frequencies,  from  this  series  in 
which  most  of  the  energy  is  concentrated. 

The  simultaneous  occurrence  of  frequency  peaks  without  clear  harmonic 
relationships  indicates  that  there  is  more  than  one  possible  mode  of 
instability  occurring  for  any  particular  flame.  There  are  several 
mechanisms  and  modes  of  combustion  instability  which  have  been  proposed  in 
the  past,  and  further  work  is  now  in  progress  to  determine  the  most 


-  17 


Important  mechanisms  in  the  present  experiments.  It  should  be  noted  that 
the  Initial  conditions  of  the  jet  flame  are  essentially  laminar  and  there 
were  no  indications,  in  the  cold  jets,  of  the  peak  frequency  values  found 
in  the  flames.  Furthermore  tests  showed  that  the  frequency  spectra  were 
insensitive  to  -the  application  of  an  external  sound  field;  this  indicates 
that  the  initial  instabilities  result  from  a  powerful  mechanism  which 
cannot  easily  be  interfered  with  by  the  imposition  of  acoustic  waves 
(due  to  fans,  instruments  etc.). 

Figure  16  shows  the  peak  frequencies  found  in  Jet  Flame  1.  The 
most  energetic  peak  was  found  to  change  to  lower  frequencies  by  discrete 
•jumps'  with  increasing  distance  downstream.  Prior  to  the  laser  schlieren 
experiments,  and  the  insight  provided  by  their  spectra,  a  series  of  schlieren 
movies  of  Jet  Flame  1  were  analysed  for  eddy  passing  frequencies  along  a 
line  above  the  nozzle  lip.  These  data  are  included  in  Figure  16  and 
the  agreement  is  seen  to  be  remarkably  good.  Thus  the  vortices  clearly 
seen  on  the  inner  side  of  the  mixing  layer  are  responsible  for  the  spectra 
peaks.  Analysis  of  cine  films  requires  subjective  judgement,  at  times,  and 
there  is  inevitably  a  'smearing  out'  of  the  coalescing,  or  other  processes, 
responsible  for  the  transfer  of  energy  between  discrete  frequencies. 

Beyond  x  =  10D,  the  dominant  frequency  is  near  30  Hz,  and  this  is  considered 
to  correspond  with  the  encroachment  of  the  large,  slow  moving,  outer 
vortices,  across  most  of  the  width  of  the  mixing  layer. 

COMPUTER  MODEL 

An  attempt  is  being  made  to  construct  a  mathematical  model  of  the 
unstable  shear  layer  in  an  axisymmetric  flame  jet.  At  present, a  vortex 
approach  is  being  used  to  model  the  large  coherent  structures  observed 
in  jets.  This  time  •dependent  model  will  be  focussed  on  the  initial  region 
of  the  jet  where  the  laminar  flow  develops  instabilities  which  eventually 
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give  rise  to  the  fully  turbulent  flow  downstream. 

The  aim  of  the  model  is  to  determine  which  parameters  are  important 
in  the  growth  of  these  eddies  and  how  and  why  the  burning  flow  differs 
from  the  corresponding  cold  flow. 

The  model  will  attempt  to  follow  the  motion  of  the  fluid  by  following 
individual  fluid  elements  in  the  flow  which  will  constitute  a  Lagrangian 
grid.  These  elements  will  be  carefully  selected  so  that  attention  can 
be  concentrated  on  the  region  in  which  mixing, heat  release  and  temperature 
fluctuations  occur,  i.e.  the  shear  layer. 

The  basis  for  the  mathematical  model  has  been  developed  from  the 
concept  whereby  a  vortex  sheet  can  be  described  by  discrete  vortices. 

This  potential  method  approach  has  been  tried  with  varying  degrees  of 
success  by  several  workers.  Although  the  flow  fields  predicted  by  these 
methods  show  many  observable  features  of  the  flow,  e.g.  roll-up  into  eddies 
and  coalescence,  they  fail  to  show  accurately  some  of  the  detail  of  the 
flow,  e.g.  the  distorted  interface;  this  is  especially  important  for 
flows  involving  combustion.  Due  to  the  sweeping  nature  of  the  approxima- 
tion  this  is  hardly  surprising.  At  present  the  model  being  used  incorporates 
some  of  the  simplicity  of  previous  vortex  models,  while  simulating  more 
of  the  'reality'  of  the  actual  flow. 

Initially  the  cold  flow  has  been  investigated,  in  which  a  cylindrical 
vortex  sheet  in  an  unbounded  domain  becomes  unstable.  To  make  the  calcula¬ 
tions,  a  periodic  disturbance  has  initially  been  given  to  the  sheet  and 
the  development  of  the  disturbance  followed.  As  can  be  seen  in  Figure  16 
the  initially  well-defined  grid  becomes  distorted  with  increasing  time 
due  to  the  build  up  of  numerical  errors.  In  order  to  overcome  this 
problem  it  may  be  necessary  to  re-define  the  grid  at  certain  stages  in 
the  flow  calculation. 

Effects  of  combustion  on  the  flow  will  gradually  be  incorporated  into 
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the  model  so  that  the  effects  of  heat  release,  expansion  of  the  gases, 
and  changes  In  viscosity  can  be  studied.  Initially  a  fast  chemistry 
reaction  will  be  assumed  to  see  if  this  gives  rise  to  the  significant 
differences  observed  in  the  burning  flow.  At  a  later  date  some  chemical 
kinetics  may  be  incorporated  into  the  model  to  see  if  this  plays  a  signifi¬ 
cant  role  in  determining  the  flow  characteristics  and  the  observed  com¬ 
bustion  instability. 

CONCLUDING  REMARKS 

Accurate  measurement  techniques,  with  good  resolution,  carefully 
designed  flows,  with  known  and  controlled  boundary  conditions,  and 
systematic  variation  of  input  variables,  have  been  used  to  gain  fundamental 
data  on,  and  insight  into,  the  structure  of  turbulent  flames.  A  double¬ 
vortex  structure  is  found:  with  outer  vortices,  which  may  originate  from 
Kelvin-Helmholtz  or  Rayleigh  type  instabilities,  and  inner  vortices  which 
are  influenced  strongly  by  one  or  more  combustion  instability  mechanisms . 
Although  most  data  so  far  have  been  derived  for  the  first  20D  near  the 
burner  nozzle,  there  are  clear  indications,  from  cine  films,  that  the 
coherent  structures  which  originate  in  this  region  can  still  be  traced 
near  the  end  of  the  visible  flames  (40D  to  100D  depending  upon  the  values 
of  $  and  Re) .  Attention  is  currently  focussed  on  the  application  of 
multiprobe  and  conditional  sampling  techniques  to  measure  the  detailed 
structures  of  the  observed  interacting  instability  waves  and  vortices. 
Because  of  the  complex  but  apparently  deterministic  nature  of  most  of 
the  flow,  more  effort  is  now  being  applied  to  time  dependent  computer 
models,  which  attempt  to  include  the  various  observed  physical  phenomena. 
The  insight  and  data  provided  by  this  study  should  enable  eventual  improve¬ 
ment  in  the  accuracy  of  predictive  techniques  for  practical  combustion 


systems . 


20  - 


REFERENCES 

1.  Chigier,  N.  A.  and  Yule,  A.  J.  "Investigations  of  turbulent 

combustion".  Interim  Scientific  Report,  Grant  AFOSR-77-3414, 
Sept.  30,  1978. 

2.  Chigier,  N.  A.  and  Yule,  A.  J.  "Large  scale  structure  in  a  turbulent 

jet  diffusion  flame  -  I",  1st  Annual  Scientific  Report, 

ONR  (Project  SQUID)  Subcontract  No.  8960-30,  April  1979. 

3.  Chigier,  N.  A.  and  Yule,  A.  J.  "The  physical  structure  of  turbulent 

flames",  AIAA  Paper  79-0217,  17th  Aerospace  Sciences  Meeting, 
New  Orleans,  Jan.  15-17,  1979. 

4.  Wohl,  K.,  Gazley,  C.  and  Kapp,  N.  "Diffusion  Flames",  3rd  Symposium 

on  Combustion  and  Flame  and  Combustion  Phenomena,  pp.  288-300, 
1949. 

5.  Ballal,  D.  R.  and  Lefebvre,  A.  H.  "The  structure  and  propagation 

of  turbulent  flames”,  Proc.  R.  Soc.  London,  A,  344,  pp.  217-235, 
1975. 

6.  Yule,  A.  J.  "Large-scale  structure  in  the  mixing  layer  of  a  round 

jet",  J.  Fluid  Mech.  89,  pp.  413-432,  1978. 

7.  Michalke,  A.  "Vortex  formation  in  a  free  boundary  layer  according 

to  stability  theory",  J.  Fluid  Mech.,  22,  2,  pp.  371-383,  1965. 

8.  Whiffen,  M.  C.,  Lau,  J.  C.  and  Smith,  D.  M.  "Design  of  laser 

velocimeter  experiments  for  turbulence  measurements". 

Laser  Velocimetry  and  Particle  Sizing,  ed.  Thompson,  H.  D. 
and  Stevenson,  W.  J.,  pp.  197-207,  1978. 


FIGURE  4.  Visual  appearance  of  Jet  Flame  1  between  10D  and  15D 


FIGURE  5.  Measurements  of  mean  velocity  in  Cold  Jet  1  by  LDA. 


FIGURE  6.  Measurements  of  mean  velocity  in  Jet  Flame  1  by  LDA.  (Re  =  10  ,  $  -  io.4) 


turbulence  intensity  in  Cold  Jet 


ionization  probe  in  Jet  Flame 
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Comparison  of  mean  velocity  and  mean  ion  current  profiles  in  Jet  Flame 


Power  Spectral  Density  of  Laser- Schlieren  Signal  (Arbitrary  Units) 
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FIGURE  13.  Spectra  of  laser  schlieren  signals 

in  vertical  line  above  nozzle  lip  in  Jet  Flame  1'. 
(Log/Log  plot) 
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Frequencies  of  peaks  in 


FIGURE  16.  Peak  frequencies  in  Laser  Schlieren  spectra  compared  with  eddy  passing  frequencies 
measured  independently,  from  high  speed  cine  films. 
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